Bimetric structure formation: non- Gaussian predictions 



Joao Magueijo^, Johannes Noller^ and Federico Piazza^ 

^Theoretical Physics, Blackett Laboratory, Imperial College, London, SW7 2BZ, UK 
^Perimeter Institute for Theoretical Physics Waterloo, Ontario, N2L 2Y5, Canada 

Abstract 

The minimal bimetric theory employing a disformal transformation between mat- 
ter and gravity metrics is known to produce exactly scale-invariant fluctuations. It 
has a purely equilateral non-Gaussian signal, with an amplitude smaller than that of 
DBI inflation (with opposite sign) but larger than standard inflation. We consider 
non-minimal bimetric models, where the coupling B appearing in the disformal trans- 
formation §1^1, = g^j^y — Bdn(j)di,(f) can run with cj). For power-law B{(j)) these models 
predict tilted spectra. For each value of the spectral index, a distinctive distortion to 
the equilateral property can be found. The constraint between this distortion and the 
spectral index can be seen as a "consistency relation" for non-minimal bimetric models. 



1 Introduction 



As more and more data pours into cosmology (e.g. [T]) the pressure is on theorists and model- 
builders to predict signatures that would unmistakably falsify cosmological models. No longer 
is it good enough to face the "zeroth-order challenge": that of producing scale-invariant 
scalar fluctuations with the correct amplitude. Deviations from exact scale-invariance will be 
detected without controversy (if they do exist) in the near future. Fields tantalizingly beyond 
our reach - such as gravitational waves (tensor modes) or primordial non-Gaussianity [2] - 
will hopefully become tangible over the next decade. While inflation [3] has dominated the 
theoretical cosmology scene, interest has never floundered on alternatives, such as pre-big 
bang cosmology [U |5], ekpyrotic/cyclic models [HI [7] and a varying speed of light (VSL) [HI 
[HI Ho]. In this paper we take recent work on VSL [HI [121 IE] one step beyond, investigating 
departures from strict scale-invariance and non-Gaussian signatures. 

Inflationary mechanisms where a varying speed of sound plays a relevant role have 
already been explored [HI [TB^j. In [15J, it was shown that an adiabatic scale invariant 
spectrum is produced even if the expansion - albeit still inflationary - is far from exponential 
(the equation of state can be as far as w ~ — 3/4 from de-Sitter), provided the speed of sound 
varies appropriately. A class of contracting ( "ekpyrotic" ) cosmologies where this mechanism 
can be applied was also found. Non-gaussianities were calculated in the limit of strict scale 
invariance (n^ = 1), they can be large in both the expanding and the contracting cases. The 
superluminal phase Cs > 1 that we are considering here allows to consider expansions that 
are not even inflationary, as long as the condition oc Cg is satisfied (here H is the Hubble 
parameter) jllj . 

Perhaps the most elegant formulation of VSL is in the guise of disformal bimetric theories, 
for which the speed of gravity differs from the speed of light [TCI [IT]. In general this is 
achieved by constructing the Einstein-Hilbert action from an "Einstein" metric g^^ (the 
Einstein frame), whilst minimally coupling the matter fields to a "matter" metric g^^v (the 
matter frame), with: 

5'=-^ / d'^xy/^ R[g^^]+ / d'^x^/~g Cm[gi,u,^Matt] + (1.1) 

in which determines the dynamics. The two metrics are related by a non-conformal 
transformation, such as: 

(jfiu = gfiu - Bd^(j)dy(j) , (1.2) 

where (f) is the "bi-scalar" field. Here B is chosen to have dimensions of M~^, so that has 
dimensions of M. (In this paper we use metrics with signature — h -|--|-, and B is defined 
so that -B > corresponds to a speed of light larger than the speed of gravity.) In the most 
general case B is an arbitrary function of 0, but in the minimal theory it's set to a constant. 

In these theories there are two light cones at any point, one for massless matter particles, 
another for gravitons. More generally the two metrics may be seen as independent represen- 
tations of the local Lorentz group (or two non-equivalent tetrads [H]), one valid for gravitons 
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and the other for matter. Thus, different Lorentz transformations must be used to trans- 
form among measurements made with matter and gravity (or equivalently, with clocks and 
rods operated by matter or gravitational phenomena). For this reason causality paradoxes 
can be skirted [l3l[T8], in contrast to straightforward tachyonic matter |^. This argument 
makes the bimetric construction important in interpreting superluminal structure formation 
models. 

A number of dynamics Sfp for bimetric theories have been considered. It was pointed 
out in [13] that a Klein-Gordon equation for in the matter frame translates into DBI 
dynamics in the gravity frame. Its corresponding Lagrangian, however, is not the Klein- 
Gordon Lagrangian in the matter frame, but simply a cosmological constant. (It was first 
noted in [T^ that for bimetric theories a Klein-Gordon action in the matter frame doesn't 
translate into a Klein-Gordon equation in that frame). Thus, the simplest bi-scalar dynamics 
is generated by 



and Am < leads to a speed of light larger than the speed of gravity. If we require the field 
to have Klein-Gordon dynamics in the Einstein frame at low energies (when matter and 
gravity frames coincide), we should consider additionally: 



i.e. a positive cosmological constant in the Einstein frame balanced by a negative one in 
the matter frame, both with magnitude tuned to 1/(25). This action maps into the DBI 
action |20l [21] in the gravity frame with DBI coupling / = ~B, as explained in [13]. For a 
choice of sign where the speed of light in the gravity frame is larger than one (/ = —B < 0), 
this is sometimes labeled "anti-DBI" (although one should note that "flipping" the sign 
of / means that this setup cannot be interpreted as portraying a relativistic probe brane 
embedded in a five dimensional bulk, as usual for DBI; for an earlier study of anti-DBI 
theories see also [22]). Combined with a mass potential in the gravity frame it leads to 
scaling solutions and scale-invariant fluctuations [231 [H] , without the need for accelerated 
expansion or a contracting pre-Big-Bang phase. The investigation of non-Gaussian features 
in models where B is allowed to "run" is the main purpose of this paper. 

In the presence of a speed of sound Cs 7^ 1 for adiabatic perturbations, the three-point 
function contains terms proportional to the power spectrum squared and terms which are 
further multiplied by a factor c^^ [211 [23 US]- In the sub-luminal case < 1, the "c^^" 
terms dominate. This is what enhances non-gaussianities in DBI inflation and makes the 
three-point function scale dependent [261 EZ] in the case of a varying speed of sound |15j 
(the combination H'^/cg is set to be constant by the scale invariance of the power spectrum 
and terms that appear with different powers of will therefore run with the scale). In the 
opposite limit, the one of infinite speed of sound that we are considering here, the "c~^" terms 
are suppressed and the remaining terms inherit the scale invariance from the power spectrum 




(1.3) 




(1.4) 
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of the two point function. The dimensionless quantity f^i is of order 1 and has opposite 
sign to DBI inflation, i.e. /tvl ~ 1 > with the WMAP sign convention. While such a 
small signal is surely observationally challenging, our results are also appealing because of 
a consistency relation between the three- and two- point functions. In fact, in the Cg ^ 1 
limit the three point function [eq. ( |4.15[ ) below] becomes independent of the background 
parameters (such as the equation of state w) and only mildly depends on the tilt — 1 of 
the power spectrum. 

The structure of this paper is as follows. In Section [2] we review and extend results 
of cosmological perturbation theory needed for the calculations in this paper. Then, in 
Section |3] we explain how scale invariance may be achieved in these models and derive 
the associated non-Gaussian features. The non-minimal model is spelled out in Section |4| 
with the basic "Gaussian" predictions presented as well as its non-Gaussian properties. 
Throughout the paper we refer to two appendices, where we explain the more technical 
aspects of the calculation. Finally in a concluding section we examine our results from a 
wider perspective. 



2 Cosmological Perturbations 



Projecting (1.4) onto the Einstein frame leads to the (anti)-DBI action, which belongs to 
the general class of k-essence models [28] . Cosmological perturbations have been extensively 
studied for these models. Here we review the main results, extending them wherever needed. 
The starting point is an action of the form: 



S 



-9 



R 



(2.1) 



where the pressure P is a general function of the scalar field and the kinetic term X 
— ^g^'^d^(f>du(f>. The energy density reads 



p = 2XP^x - P : 



while the speed of sound is given by 



P 



X 



P,x 



Px + 2XPxx 



(2.2) 



(2.3) 



In a FRW Universe of scale factor a{t) and Hubble rate H{t) = a/ a we define the slow-roll 
parameters as follows: 



H 
IP' 



V = 



Vs 



(2.4) 



In the Us = 1 scale-invariant case, non-Gaussianity has been calculated in [U]. Here we 
generalize to the case of an arbitrary - albeit small - tilt and negligible running [r] ^ 



0,rjs ^ 0). The calculation of the three-point function also necessitates defining two further 
parameters derived from P{X, (p) [25] 



A 



XPx + 2X^Pxx 



X P,xx + 2"^ P,xxx 



(2.5) 
(2.6) 



At quadratic order, the action for the curvature perturbation ( for general speed of sound 
models is given by [29] 



^2 



^ / d'xdr 



dC 
dr 



(2.7) 



where r is conformal time, and z is defined as usual by z = aylejcs- In [15] it was shown 
that it is convenient to work in terms of the "sound-horizon" time dy = c^dr instead of r. 
Explicitly, when 77 = 77^ = 0, 

1/ = T^-^^TTTTT • (2.8) 



It is useful to write the behavior in |/-time of some relevant quantities: 



-2/)' 



E 



-y)' 



The quadratic action then takes the form 

S=^ j d^xdy q 

where ' = d/dy, and 

q = ^/c^z = 



C - (VC)^ 



aV2e 



(2.9) 



(2.10) 



(2.11) 



Upon quantization the perturbations are expressed through creators and annihilators as 
follows, 

ay,k)=Uk{y)a{k) + ul{y)a\-k). (2.12) 

However, in order to make the correct choice of vacuum, it is useful to refer to the canonically- 
normalized scalar variable v = Mp\q(. Then the equations of motion for the Fourier modes 
are given by 



Vk = 0. 



(2.13) 



It is well-known that this results in a scale-invariant spectrum if q" /q = 2/y'^. More generally, 
we have 

q" If, l\ ..... 



and the solution for Vk{y) corresponding to the Bunch Davis vacuum is 



Vk{y) 



-y 



Hi'\-ky), 



(2.15) 



where hI^^ are Hankel functions of the first kind. The relation between ris, ^ and e and is 

2e + e, 



3 -2z/ 



e- 1 



Again, in the limit when v = 3/2, 

Vk{y) = -- 



1 



ky 



-iky 



[u = 3/2) 



(2.16) 



(2.17) 



and we recover the scale invariant spectrum. Going back to the modes defined in (2.12) we 
have 



Uk{y) 



.1/2 



aMpiV2e 



Vk{y) 



.1/2 



aMpi23/2 



yHi'\-ky). 



(2.18) 



It is useful to adopt an approximate expression for the Hankel functions. By expanding at 
\ky\ 1 we obtain 



Hi'\-ky) = -t ?!£MLM:![i + ,ky + 0{kyf]e-'^y 



which gives the following approximate expression for Uk, 



Uk{y) 



2Mpiy/^J^ \ 2 J 



3/2-1/ 



;^1 + iky)e 



-iky 



In order to obtain (2.20), eq. (2.8) has been used, together with V{i> ^ 3/2) 



(2.19) 



(2.20) 
/2|3 As 



expected, Uk{y) — const, in the y ^ Q limit. To check this explicitly we use ( 2.9[ ) and note 
that 

— - <-yr-'" . 



.1/2 



The derivative of Uk{y) with respect to y is also easily obtained: 



u'kiy) 



.Hie + Es-l] 



-ky 



k'^y e 



-iky 



2Mpiv/^;Fe 

Finally, the expression for the C Power Spectrum reads 

1 (e, + e-l)2 22--3 



C 



27r 



^k'\Ck\ 



2(27r)2e c.M^j ' 



(2.21) 



(2.22) 



(2.23) 



where the bar symbol means that the corresponding quantity has to be evaluated, for each 
mode fc, at sound horizon exit, i.e., when y = k~^. 

iln fact, T{iy) = 0f/2 [1 + 0.036(i/ - 3/2) + ...]. 
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3 The Scale Invariant limit 



For the bimetric theories discussed in the introduction, in the Einstein frame the action takes 
the form: 

p(x,0) = -ri(0)v/i-2/(0)x + ri(0) - (3.1) 

where X = —kd„(l)d^(j). Scahng solutions of this action have been studied in [301 113 [IS]- 



particular, a scale invariant spectrum of primordial perturbations is produced if 

= - 

V{<j)) = Vc 



3 


1 1 




4 " ^ 




6 V 1 


'UpJ [ 


Vo 


\MpJ 



Mp 
2e 



(3.2) 



3 1 + 
4 

3+-- 



8eMi 



(3.3) 



where we have made use of the slow-roll parameter e = —H/H^. In the large (f) limit (which 
also corresponds to the large Cs limit) we recover the form discussed in [13]: /(0) ^ —B < 0, 
V{(j)) ~ (fP'. This corresponds in general to Cg oc p. In such a strict — )■ cxd limit the 
amplitude of the three point function A can be read straightforwardly from ( |B.12 ) (we refer 
the reader to Appendices A and B for the more gruesome technical details). Comparing 
with the cubic effective action (A.l) we find that only the ^^^2 and Ac^i^oc,)'^ terms are not 



subdominant as — )■ 00. The resulting total amplitude is independent of the parameters 
{w or e) and reads 



Ac 



1 



PP - 



(3.4) 



This is precisely the equilateral shape, peaking for ki = k2 = k^, that is also obtained in the 
scaling solutions considered in [15] in the e — )■ 0, a — limit. More specifically one obtains 



A, 



e-S>0 



1 - 



1 



Ar 



+ 0(n,-l) 



(3.5) 



This shows that in the minimal bimetric model the dimensionless quantity Jml is of order 1 
and has opposite sign to DEI inflation, i.e. fiyL ~ 1 > 0. Thus the model is quite distinct 
in this respect to standard inflation (for which f^L ~ e ~ 0.1) and DBI inflation (for which 
Inl ~ —100 is a distinct possibility.) Notice there's been some confusion [21], both among 
theorists and observers, regarding the sign of Jnl- Here we adopt the convention used by 
WMAP, where positive /nl physically corresponds to negative-skewness for the temperature 
fluctuations: we assign a negative /atl to DBI inflation, so that /atl > for the anti-DBI 
models under consideration. 
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4 Beyond the minimal model 



It could be that the parameter B appearing in the disformal transformation (1.2) is itself 
a function of 0. In this Section we show that non-minimal theories with power-law -8(0) 
lead to tilted spectra, without running. Naturally, more complicated -6(0) would lead to 
more complex spectra, so one can't say that absence of running is a general feature of these 
models. 

First we note that the speed of sound in the gravity frame is: 



K 



X 



' K^x + 2XK^xx 
whereas the density and pressure are: 

1 



2BX 



(4.1) 



P 
P 



2XKx - p 



B 



K-V 



1 

B 



+ V 



V 



(4.2) 
(4.3) 



Thus scaling solutions may be obtained for a variety of potentials, with the property that 

H 3 



-(1 + w) 



(4.4) 
(4.5) 



are constant. As explained in the Appendix A of reference jT5] the constancy of e and 
permits a simple integration into Cg = Cs{4>) and H = H{(j)). The Friedmann equations then 
give a solution for V(0) and /(0). 



V 



PI 



Ml 



PI 



1 - 



2e 



3 1 + 



3e202 



Ml 



PI 



f{<P) 



12 



3e2 



1 

2+^ 



64e2M^i 



l + O 



M 



PI 



(4.6) 
(4.7) 

(4.8) 
(4.9) 



Using B 



-/ > and in the limit Cs ^ 1 we therefore obtain 
V = Vol 



M 



PI 



(4.10) 



7 



B 



-f 



2+ 



16e2\/o \Mpi 

Following the calculation in [TT] we find that for these solutions the spectral index is: 

e, + 2e 

ns-l = — . T ■ 

es + e-1 

Whilst scale-invariance is associated with the universal law 



(4.11) 



(4.12) 



Cs (X p 



(4.13) 



and also / = 1 and a quadratic potential for all equations of state, the same doesn't happen 
if we depart from scale-invariance. Indeed exact scale invariance requires = 2e so that 
these parameters fall out of conditions for the spectral index; but this doesn't happen as 
soon as 7^ 1. Note, however, that any scaling solution has: 



(4.14) 



an expression that will be essential in evaluating non-Gaussianities. 

In Appendix B we present expressions for the Non-Gaussian amplitude A for general Cg 
profiles. Our calculations there consequently apply to both subluminal and superluminal 
cases, as long as scaling solutions with constant e and are considered. Once a solution 
for Cs is provided, these uniquely determine the Non-Gaussian signature. However, here 
we exclusively focus on the Cg — t- oo case relevant in the bimetric context. Considering 
further the relation (4.14), fixing Cs{(f)), we have that in the large limit the appropriate 



Non-Gaussian amplitude to compute is still ^5^-5-00 • 

Combining terms from the general results of Appendix |B| eqs. (B.12) and (B.13) and 
taking the small tilt (n^ 



1^1) and Cs — )■ 00 limits, we find 



A 



+ ins-l 



ris-l 



+0 



1 



(4.15) 



where the only dependence on e and appears either in the "observable" combination 
— 1 or in the subleading (9(1/Cs) terms. Upon approaching scale invariance only the first 
line inside the square brackets stays relevant, ensuring that A reduces to the equilateral 



amplitude (3.4) as required. 
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(a) -^(l,a;2,2;3)/(a;2a;3) for = 1 



(b) ~A{l,X2,X3)/{x2X3) for Us = 0.96 



Figure 1: We plot the non-Gaussian amplitude from Eq. (4.15) —A{l,X2,X3)/{x2X3) for 
Us = I (left) and Ug = 0.96 (right). 



The amplitude (4.15 ) is plotted in Figure 1 and peaks in the equilateral limit ki = k2 = 



In the local limit ki <^ A;2) ^3; on the other hand, the first line inside the square brackets of 



Eq. (4.15) goes to zero. In agreement with the consistency relation [3TI 132] we then have 

1 , 



1 



(4.16) 



The predictive power of our result lies in establishing a consistency relationship between 
ris and A. In fact, we find a distinctive Non-Gaussian signal for any given spectral index n^. 
Whilst the overall non-Gaussian amplitude A still peaks in the equilateral limit ki = k2 = k^ 
in both red- and blue-tilted cases, its shape is modified when compared with the scale 
invariant limit. Illustrating this point. Figure la shows A for the exactly scale invariant case 
Us = 1, whereas in Figure lb we plot A for a red-tilted power spectrum with = 0.96. 
Specifically we find an A^Ug = 0.96) = A{ns = 1) + A^, where is approximately one 
order of magnitude smaller than A{ns = 1). 



5 Conclusions 

Strict scale-invariance has been associated with superluminal bimetric models, where the 
speed of light is larger than the speed of gravity in the early Universe [HI |13]. Indeed 
this is a feature of the minimal bimetric model, but in this paper we showed how tilted 
spectra, red or blue, could be generated by a non-minimal bi-scalar coupling B{(f)). At first 
this might suggest we've fallen into the "theory of anything" trap, but it's not the case. A 
unique non-Gaussian shape is predicted for any value of the spectral index, with distinct 
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distortions away from the scale-invariant equilateral shape appearing for each of the tilted 
cases. These distortions can be seen as "consistency conditions" for this class of models. 
This is particularly relevant given the absence of gravitational waves for all bimetric models 
of this kind. (Note that these models solve the horizon problem for matter but not for 
gravity, so tensor modes don't start their lives inside the horizon.) 

One might wonder where the proposed running coupling B{(f)) comes from. First note that 
we don't need the full (anti-)DBI action (3.1) resulting from (1.4), unless we impose Klein- 
Gordon dynamics for (p in the Einstein frame at low energies. This may not be necessary, 
and if we relax this requirement all we need is (1.3), i.e. a negative cosmological constant 
Am in the matter frame (which, we stress, does not lead to an AdS solution). In fact, if 
we relax the low-energy requirement, doesn't even need to be related to B. If, however, 
we do insist on Klein- Gordon dynamics for (p in the Einstein frame at low energies, then 
the negative matter frame cosmological constant should be exactly balanced by a positive 
Einstein frame cosmological constant, and their common magnitude should be 1/{2B). 

A number of interesting theoretical connections can be made. In the context of emer- 
gent geometry, it's been pointed out that different emergent metrics may apply to bosons 
and fermions [53]. The fact that the vacuum energy is negative for fermions and positive 
for bosons suggests an action of the proposed form, with a speed of light larger than the 
speed of gravity (i.e. an anti-DBI action in the Einstein frame). Also these models become 
asymptotically a cuscaton [23j model, a feature that may be used to support the view that 
they are a UV-complete alternative to inflation. Finally, it is possible that this construction 
results from an entirely different set up, such as deformed dispersion relations [51]. It is 
interesting that the dispersion relations needed for scale-invariance are of the same form as 
those discussed in the context of Horava-Lifshitz theory [35]. More generally a connection 
with deformed special relativity remains to be fully explored [Ml EZ]. Absence of exact 
scale-invariance could then be a major clue into the foundations of these theories. 

While work on these theoretical ramifications is an interesting motivation, and should be 
pursued further, in this paper we focused on the phenomenology of these models. Measuring 
the shape of the three-point correlator (as opposed to a quantity as muddled as /atl) poses 
an interesting observational challenge. The fact that the matter appears coupled to the 
measurement of makes these models an interesting target for future experimental work. 
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A Appendix: The cubic action 

It is useful to report here the cubic effective action derived in [211 125] . The result is valid 
outside of the slow-roll approximation and for any time-dependent sound speed: 



S I 



Mp,2 / dM^x|-a^ 



+ 2A 



+ -(e - 2e, + 1 - ci)adCy - 2a-adOidx) 



+ 



a^e d f T] 



CC + 7^{dO{dx)d\ + -rid'Oidxf + 2/(0 



2a 



4a 



5L 



(A.i: 



where dots denote derivatives with respect to proper time t, 9 is a spatial derivative, and x 
is defined as 



(A.2) 



Meanwhile, in the last term denotes the variation of the quadratic action with respect 
to the perturbation 



5L 



(A.3) 



/(C) 



-C + —CC + — — 



V ,2 



-(90(90 + 9-^(9,9,(9,C9,c))] 



+ 



\{dQ{dx)-d-\d^^,{d,Cd,x))\ , 



(A.4) 



where 9 ^ is the inverse Laplacian. Since f^|i is proportional to the linearized equations of 



motion, it can be absorbed by a field redefinition 

C -> Cn + /(Cn 



(A.5) 



B Appendix: The three-point function 

The three point function can be calculated by following the same method of ^15j and gener- 
alizing it. The standard calculation [211 12S], at first order in perturbation theory and in 
the interaction picture, leads to 

(C(t,ki)c(t,k2)c(t,k3)) = -I fdt\[atM)at,k2)atM),H,nt{t')]) , (b.i) 

J to 



where i^int is the Hamiltonian evaluated at third order in the perturbations and is directly 
derivable from (A.l) and vacuum expectation values are evaluated w.r.t. the interacting 
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vacuum \Q). By using (2.12) and applying the commutation relations [a(k), a'''(k')] = 
(27r)^(5^(k — k'), we can calculate the three point function for each term appearing in the 
action (A.l). It is useful to follow in detail the calculation for the "CC^" piece. We have: 



(C(ki)C(k2)C(k3))(.^2 = i(27r)^5^(ki + k2 + k3)Mfcj(?/cnd)Mfc2(l/end)Mfc3(l/cnd) 

X / dy e - 3 + 3cjufc^ y — f -f + perm. + c.c. B.2 

J-oo+ie a cj ' dy dy 

The subscript "end" means that the quantity has to be evaluated at the end of "inflation" . 
We now substitute (2.20 ) and use (2.21 ) to take some time-independent combinations outside 
the integral and evaluate them ai y = yend- 

(C(k,)C(k.)C(k,))«. = ,:(2.)^^^(k.H-k..H-k3) "'f - ""'^y 



/Vend g2 
dy{t - 3 + 3c^) — (1 - ikiy)klkly'^e'^y + perm. + c.c. , (B.3) 
■oo+ie 



where we have dropped a factor of nj(l + ikjyend)^ ^^yend g^g ^j^jg i^q negligibly small i 
the limit \ky\ << 1 where the truncated Hankel function expansion (2.19) is valid. 



i nmn \Ky\ i wnere tne truncatea nauKei mnction expansion [z.iif) is vaiia. 
By using ( |2.9 ) we can finally express the time dependent quantities inside the integrals 
as power-laws in y. It is useful to report some of the basic results of [15] for integrals of this 
type. By calling 



C 



rvcnd f y ^ ^ 


—oo+ie VZ/cnd 



-zyre 



n iKy 



(B.4) 



For 7 + n > —2 the imaginary part of (B.4) is convergent as y^nd 0. In this case we can 
approximately extend the upper limit of integration to 0, which amounts to neglecting terms 
of higher order in (fc|?/end|)- We thus obtain 



ImC = -(ir|i/e„d|)-^ cos ^r(l + 7 + 



The two types of behavior that we encounter are, in particular, 



end 



^ _y_ 

1 1 



[1 - e - ^sYH'^Cs (1 - e - esYHl^^Csend Vl/cnd 



«2 



where 



02 



— 1 = 3 — 2z/ 

2e-e, 
+ e - 1 
12 



2e + e, 
+ e - : 



(B.5) 

(B.6) 
(B.7) 

(B.8) 
(B.9) 



By using the above formulas and re-expressing everything in terms of quantities calculated 
at sound horizon crossing [i.e. when, by convention, y = K~^) we finally obtain 



(C(ki)C(k2)C(k3)),/2 = (27r)^5^(ki + k2 + k: 



xMM;(e_3)cos^r(i + «2) 



1 + (1 + a2) 



K 



+ 3c^ cos + ai) 



1 + (1 + ai) 



K 



+sym. 
(B.IO) 



The three point function is conveniently expressed, after factoring out appropriate powers 
of the power spectrum, through the amplitude A, 



(C(ki)C(k2)C(k3)) = (27r)^(53(ki + k^ + \^,)P^^^^A. 



(B.ll^ 



Again, by convention, the power spectrum Pf in the above formula is calculated for the mode 



K. By following the same steps above for each of the terms in the action (A.l) we obtain 



3~2u 



1 / kik2h 



[(e - 3)X^^2(a2) + 3c^X^^2(ai)] 



A, 



1 f kik2k^ 



8c2 V 



3-2u 



[(e - 2es + l)X^(a^)2(a2) - c^^ac)^' 



A. 



A. 



1 / kik2k 



3-2u 



3-2v 



1 / kik2h 



IQcl V 2K^ 



[e^XA<^2)\ 



(B.12) 



where 



X(*^2(«) = cos ^^r(l + a) 



avr, 
~2 



^aaoM) = - cos ^r(l + a)[Y^ kj 



K 1 



l + a. 



a — 



: COS —1 (1 + a) 



fcjfc- - akik2k3 
1 — a) ■' 



1 _ 2^ ^- + 7^ '^i ^- ;^2 S- 



13 



an 



.a 



cos — r(l + a) 



X,2(a) =cos— r(l + a)(2 + a/2) 



(B.13) 



In the above, A^^ accounts for the d(dxd'^x ^'iid {d'^QidxY terms and the first term in (A.l ) 
has not been considered because it is identically null in the case of a DBI-type action. 

One should note that, in contrast to the — t- oo limit relevant in the bimetric context, 
in the subluminal case [cg < 1), it is in fact the "c^^" terms that dominate. 
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